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We aim at the construction of dark energy models without exotic matter but with a phantom-like 
equation of state (an effective phantom phase). The first model we consider is decaying vacuum 
cosmology where the fluctuations of the vacuum are taken into account. In this case, the phantom 
cosmology (with an effective, observational ui being less than — f ) emerges even for the case of a 
real dark energy with a physical equation of state parameter u larger than — f . The second proposal 
is a generalized holographic model, which is produced by the presence of an infrared cutoff. It also 
leads to an effective phantom phase, which is not a transient one as in the first model. However, we 
' show that quantum effects are able to prevent its evolution towards a Big Rip singularity. 
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lO ■ I- INTRODUCTION 

o ■ 

Recent observational constraints obtained for the dark energy equation of state (EOS) indicate that our 
Universe is probably in a superaccelerating phase (dubbed as "phantom cosmology"), i.e. H > (see 
Oh. and references therein). When dark energy is modeled as a perfect fluid, it can be easily found that in order 

to drive a phantom cosmology, the dark energy EOS should satisfy = Pd/Pd < —1 (phantom) 0|. The 
classical phantom has already many similarities with a quantum field However, it violates important 
energy conditions and all attempts to consider it as a (quantum) field theory show that it suffers from 
instabilities both from a quantum mechanical and from a gravitational viewpoint 0. Thermodynamics of 
this theory are also quite unusual ^ the corresponding entropy must be negative or either negative 
' temperatures must be introduced in order to obtain a positive entropy. The typical final state of a phantom 

universe is a Big Rip singularity ^ which occurs within a finite time interval. (Note that quantum effects 
coming from matter may slow up or even prevent the Big Rip singularity which points towards a 

transient character of the phantom era). 

In the light of such problems, we feel that it would be very interesting to construct a new scenario where 
an (effective) phantom cosmology could be described through more consistent dark energy models with 
ujd > —1. Recently, some proposals along this line have already appeared. An incomplete list of them 
includes the consideration of modified gravity |p or generalized gravity |lOj |. taking into account quantum 
effects non-linear gravity-matter couplings . the interaction between the photon and the axion 13], 
two-scalar dark energy models *he braneworld approach (for a review, see j^, however, FRW brane 
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cosmology may not be predictive jig); and several others. 

In this paper we continue the construction of more consistent dark energy models, with an effective phan- 
tom phase but without introducing the phantom field, what is realized by using a number of considerations 
from holography and vacuum fluctuations. In the next section, a dark energy model which includes vacuum 
fluctuations effects (a decaying vacuum cosmology) is considered. With reasonable assumptions about the 
vacuum decaying law, the FRW equations are modified. Then, it is shown that the effective dark energy 
EOS parameter, which is constrained by astrophysical observations, is less than —1. However, even if the 
dark energy EOS parameter is actuall bigger than —1, the modified FRW expansion, or in other words, a 
decaying vacuum can trick us into thinking about some exotic matter with w less than —1. Hence, such 
dark energy model does not need the exotic phantom matter, its effective phantom phase is transient and 
it does not contain a Big Rip singularity. In section 3 another holographic dark energy model which uses 
the existence of an infrared cutoff, identified with a particle or future horizon, is constructed. It generalizes 
some previous attempts in this direction and leads to modified FRW dynamics. Such model also contains 
an effective phantom phase (without introduction of a phantom field), which leads to Big Rip type singu- 
larity. However, taking into account matter quantum effects may again prevent (or moderate) the Big Rip 
emergence, in the same way as in Refs. ^^]. A summary and outlook is given in the Discussion section. 



When considering the energy of the vacuum, there are two conceptually different contributions: the energy 
of the vacuum state and the energy corresponding to the fluctuations of the vacuum. More precisely, the 
first is the eigenvalue of the Hamiltonian acting on the vacuum state: E — {0\H\0), while the second is the 
dispersion of the Hamiltonian: Ea = (0|i?^|0) - {{0\H\0})'^. 

In constructing the dark energy as a scalar field ("quintessence") only the energy of the vacuum 
state is actually considered, and the vacuum fluctuations are neglected. Note that, while the potential 
energy of the scalar field corresponds to the energy of the vacuum state, the kinetic energy density does 
not correspond to the energy of the vacuum fluctuations. The kinetic energy corresponds to the energy of 
the excited states. Actually, ignoring vacuum fluctuations may not be a well justified assumption. It was 
ignored for quintessence models, just for simplicity. 

Thus, confronted with the difficulties of building a phantom universe from quintessence alone, it is natural 
to consider the sector which was ignored before. In this paper, we will show that the ignored sector can 
really help to resolve the problems above. The key point is that, when considering those two contributions 
to the vacuum energy in a cosmological setup, their energy density changes with the cosmological expansion 
in different ways. 

First, if dark energy is modeled as a scalar field 4>, interacting only with itself and with gravity ^|, then 
it behaves just like an ordinary perfect fluid satisfying the continuity equation 



II. DARK ENERGY AND VACUUM FLUCTUATIONS 



Pd + 'iH{pd + Pd) = 0, 



(1) 
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Pd 




+ 



(2) 



and 



Pd 




(3) 



Thus 



UJd = 



> -I. 



(4) 
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When comparing dark energy models with observations, it is often assumed that ujd = const, (this is reason- 
able since allowing EOS to vary will increase the number of parameters to constrain). Then from equation 
(P) it follows that 

PdOc(l + z)3(i+-''). 

On the other hand, to find the energy corresponding to the fluctuations of the vacuum, we must have a 
regularization scheme in order to handle the divergent expression [3. Since the regularization scheme must 
be Lorentz invariant (as is the case, e.g. of dimensional regularization in curved spacetime), it turns out 
that the final result must be of the form Tf^^^, — pAgfif even if is time- varying due to the time variation 
of the IR cutoff 19]. Note that mathematically, ptvg^p is exactly the energy- momentum of a perfect fluid 
with p = —p. Then from the energy-momentum conservation law, we find that the vacuum fluctuations will 
interact with matter [20 ll2lll2^ 

Pm + 3Hp,n = -pA, (6) 

where pm is the energy density of nonrelativistic matter and pA is the energy density of vacuum fluctuations. 
Following ^27] , we will call a cosmological model which includes the effects of vacuum fluctuations a "decaying 
vacuum cosmology" since, as can be seen from ©, during the expansion of the Universe, owing to the fact 
that Pa is decreasing, new matter is created from the vacuum. 

It is clear that Eq. ^ and the Friedmann equation are not enough to describe a decaying vacuum 
cosmology. We should also specify a vacuum decaying law. There are lots of papers on various proposals 
concerning how p^ will change with cosmological time (see 0| and references therein). However, recently 
it has been proposed that from 10) and a simple assumption about the form of the modified matter 
expansion rate, the vacuum decaying law can be constrained to a two-parameter family. Thus, we can think 
about the observational consequences of a decaying vacuum cosmology, even if we do not understand the 
physics underlying the decaying law. More precisely, the main (and rather natural) assumption is that the 
matter expansion rate be of the form 

Pm = Pmoa^^^^ ■ (7) 

where p„iO is the present value of p„j. Actually, this assumption is valid in all the existing models of a 
decaying vacuum cosmology. An immediate simple observation is that one must have e < 1. Otherwise, 
the Universe would expand with an acceleration in the matter dominated era, which is excluded by the 
observations of SNe la, that our Universe expanded with deceleration before redshift z ~ 0.5 1]. Actually, 
it is expected that e <C 1 since, so far, there has been no observational evidence about an anomalous dark 
matter expansion rate. 

From (O and {Tj), we can find that pA is given by 

PA = PAO + -z a , 8 

where pAo is an integration constant. In 21] it was shown that all the existing vacuum decaying laws can 
be described by a suitable choice of e and pA- In our current framework, the term pA can be absorbed into 
the dark energy density, so in the following discussion, we will simply assume pA — 0. The FRW equation 
describing a Universe consisting of cold dark matter, dark energy and vacuum fluctuations is 

H^ = ^{Prn+Pd + PA)- (9) 

It is important to note that due to the modified dark matter expansion rate, the physical dark energy EOS 
Ud is no longer the one directly seen in supernova observations. To compare our model with observations, 
we can adopt the framework of effective dark energy EOS 23] that can unify quintessence, modified gravity, 
and decaying vacuum cosmology into one and a single framework. The effective dark energy EOS is defined 
as 111 

, IdlnSH^iz) 
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where 5H^ = — Pmo(l + zY / i^Mp) characterizes any contribution to the cosmic expansion in addition 
to the standard cold dark matter. It is oJeff that is the real quantity which is constrained in cosmological 
observations. The explicit form of 5H^ in the above model is 

from which we can find the effective dark energy EOS: 



From this one can see that, ignoring the decay of the vacuum, i.e. with e = 0, then oJeff = ^d- 

FIG. n shows the evolution of We// up to redshift 2 for e = 0.01, = 0.9 from top to bottom. It is easy 
to see that LOeff — for small z and negative z. As mentioned above, in current works constraining i^eff, 
the most reliable assumption is that uJeff is a constant 0. Thus, what is actually constrained in those 
works is the average value of the effective EOS We// = / ^ef f{o)uJef f{a)da/ J VLf.f f{a)da p^ . which from 
FIG. n is smaller than —1 if averaged up to redshift 2 (the current upper bound on observable supernova) 
for roughly e > 0.05. Thus, we have constructed a model which can drive a phantom cosmology with a 
physical EOS larger than —1. In this picture, uieff < — 1 is not the result of exotic behavior of the dark 
energy, but the modified expansion rate of the CDM, i.e. a decaying vacuum can trick us into thinking that 
UJeff < —1- Under the assumption of a constant EOS, the current constraint on uieff from SNe la data is 
-4.36 <0Jeff < 0.8 at 95% C.L. Q. Then e < 0.1 is obviously consistent with current constraints on dark 
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energy EOS. Moreover, it is easy to see from FIG. that ujeff will be more and more negative when the 
redshift gets larger. Thus, if more high redshift SNe la data become available in the future, the interaction 
between vacuum and nonrelativistic matter predicts that we will get a more negative value for i^eff- 

It is also interesting to note that in the current model, the arguments leading to a Big Rip of the Universe 
in a phantom cosmology no longer apply. While the current expansion rate is superaccelerating, the 
Universe will be driven by only the healthy ljx > —1 dark energy in the future. The Universe is now 
superaccelerating just because non-relativistic matter is of the same order of magnitude as dark energy, 
and thus the effects of the nonstandard expansion rate of nonrelativistic matter cannot be ignored. This 
conclusion is supported by recent analysis on the construction of a sensible quantum gravity theory in 
cosmological spacetime |23. Actually, it is easy to see that the trouble with de Sitter space discussed in 
[25j| will be avoided in a superaccelerating cosmology: the total amount of Hawking radiation due to the 
cosmological horizon is infinite and thus it will thcrmalize the observer. Furthermore, since the horizon is 
decreasing in a superaccelerating cosmology, black holes will dominate the spectrum of thermal fluctuations 
in late time [2^. An observer may be swallowed by a horizon sized black hole before he/she can reach the Big 
Rip singularity. Thus, it is doubtful that the appearance of Big Rip singularity makes sense in the quantum 
gravity theory. This is also supported by a direct account of the quantum gravity effects moderating or 
preventing the Big Rip 0, Q . 

In the above discussion, the simplest model of dark energy characterized by a constant EOS is considered 
with the assumption that the vacuum decays only to dark matter. Let us now discuss a more general case. 
Let the vacuum fluctuation energy decay into both dark matter and dark energy. Specifically, we will consider 
a dark energy with EOS proposed in (see, also [EH^I), 

Pd = -Pd-Ap'^, (13) 

where A and a are constant. This can be viewed as characterizing the minimum deviation of the dark energy 
EOS from that of a cosmological constant. When considering the effect of vacuum fluctuations, from the 
total energy-momentum conservation, one gets 

Pd + Pni+iH{pd+Pd + Pyn) ^ -PA- (14) 

Now one cannot proceed further without specifying a vacuum decaying law. Currently, the following law 
can be proposed 

PA = eMlH^. (15) 

This law can be deduced from the results of |18|. in which the vacuum fluctuations of a spatially finite 
quantum system, e.g. free scalar field, have been evaluated. The most interesting conclusion in this compu- 
tation is that the energy corresponding to the vacuum fluctuations is proportional to the area of the spatial 
boundary of the system. In a cosmological setting, it is reasonable to view the Hubble horizon as the spatial 
boundary of the quantum theory describing matter in the Universe. 

From equation (|15|) and the Friedmann equation, it can be shown that Eq. H14|l can be written as 

Pd + Prn + (3 - e)H{pd +Pd + Prn) = 0, (16) 

from which we find that 

P<i-p<io[l-(3-e)i(l-a)ln(l + z)]i/(i-"), = p™o(l + ^)'"% (17) 

where A = For e = 0, i.e. in the absence of vacuum decaying effects, (|17|l reduces to the expression 

found in 8, 26]. 

Now, one can find the effective dark energy EOS by exactly the same procedure as above: 

(1 + - (1 + ^? - ^Vol"^" (1 - (3 - e)i(l - a) ln(l + z))"/(i-") 
'^"^ 3^(1 + z)3-^ -(1 + ^)3 -1-^(1- (3 -e)i(l- a) ln(l + z))i/(i-a) ' 
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FIG. 2: Evolution of the physical (thin line) and effective (thick line) dark energy EOS as given by equations H19|l 
and TM . 

When e = 0, this reduces to the physical equation of state, 

uJd = — = -l-A{l- 3A(1 - a) ln(l + z))-\ (19) 
Pd 

FIG. 121 shows the evolution of ujeff and uj^ up to redshift 2 for yl = —1/3, a = 1/2, e = 0.1. Hence, 
it is obvious that while uj for the physical dark energy EOS is larger than —1, with account of vacuum 
fluctuations, for the dark energy EOS is smaller than —1. Thus, we have constructed a vacuum fluctuation 
dark energy model where the phantom phase is transient and there is no need to introduce exotic matter (a 
phantom) to comply with the observational data. 

III. HOLOGRAPHIC DARK ENERGY AND THE PHANTOM ERA 

The fluctuations of the vacuum energy are caused by quantum effects. It is expected, that the quantum 
corrections to the theories coupled with gravity should be constrained by the holography. The holography 
requires an infrared cutoff L\, which has dimensions of length and could be dual to the ultraviolet cutoff. 
Hence, the fluctuations of the vacuum energy are naturally proportional to the inverse square of the infrared 
cutoff. This points to the connection of dark energy based on fluctuations of the vacuum energy with the 
holographic dark energy. Our purpose here will be to construct the generalization of the holographic dark 
energy model which naturally contains the effective phantom phase. 

Let us start from the holographic dark energy model which is a generalization of the model in [23| (see 
also [2^ [2^ 1131 )■ Denote the infrared cutoff by La, which has a dimension of length. If the holographic 
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dark energy pA is given by, 

3c' 

with a numerical constant c, the first FRW equation (H'^/k^) — p\ can be written as 



PA = , (20) 

ft, i^A 



"'11- 

Here it is assumed that c is positive to assure the expansion of the universe. The particle horizon Lp and 
future horizon L f are defined by 

dt ^ f°° dt 

Lp = a — , Lf = a — . 22 
Jo a Jt a 

For the FRW metric with the flat spacial part: 

ds^ = -dt^ + a{tf ^ {dx^y . (23) 

4=1,2,3 

Identifying L\ with Lp or Lj, one obtains the following equation: 

dt i^aH^ a ^ ^ 

Here, the plus (resp. minus) sign corresponds to the particle (resp. future) horizon. Assuming 

a = aot^° , (25) 

it follows that 

ho = lTc. (26) 

Then, in the case La — Lf, the universe is accelerating {ho > or w = — 1 + 2/3ho < —1/3). When c > 1 
in the case L\ = Lp, ho becomes negative and the universe is shrinking. If the theory is invariant under the 
change of the direction of time, one may change t with —t. Furthermore by properly shifting the origin of 
time, we obtain, instead of \2b\i . 

a = ao{ts-tt° . (27) 

This tells us that there will be a Big Rip singularity a,t t — ts- Since the direction of time is changed, the 
past horizon becomes a future like one: 

Lp^Lf^aj - = al — . (28) 



a 



Note that if La is chosen as a future horizon in H22() . the deSitter space 

a = ancT = (29) 

can be a solution. Since Lf is now given by i/ = /, we find that the holographic dark energy H20|l is given 

by PA — -^rp- Then when c = 1, the first FRW equation -^Ll^ — Pa is identically satisfied. If c ^ 1, the 
deSitter solution is not a solution. If we choose La to be the past horizon, the deSitter solution docs not 
exist, either, since the past horizon Lp ((^ is not a constant: Lp = (1 — e")//. 
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In [SJ, it has been argued that, when we choose La to be the particle horizon Lp, the obtained energy 
density gives a natural value consistent with the observations but the parameter w of the equation of the 
state vanishes since Lp should behave as a(i)'^/^. That seems to contradict with the observed value. In the 
present paper, we wish to generalize the holographic model to have, hopefully, more freedom and to be more 
successful in cosmology. As we will see in this section, there may be many possiblilities for the choice of 
the infrared cutoff L^ which are more consistent with observations (giving correct cosmology without above 
problems). 

In general, L\ could be a combination (a function) of both, Lp, Lf. Furthermore, if the span of life of 
the universe is finite, the span tg can be an infrared cutoff. If the span of life of the universe is finite, the 
definition of the future horizon L f (|22() is not well-posed, since t cannot go to infinity. Then, we may redefine 
the future horizon as in 

dt f°° da 



Lf^Lf^aj -^aj —. (30) 







Since there can be many choices for the infrared cutoff, one may assume La is the function of Lp, L^, and 
is, as long as these quantities are finite: 



La = La (^Lp, L/, . (31) 



As an example, we consider the model 



La 



2t. 



c r / , , ^ X 2 



(32) 



which leads to the solution: 



In fact, one finds 



I da f " dt t , , 

L, + L,=al j^-a - ^ (34) 



-f-f" Jo 
and therefore 

c 1 /I 1 



La 2 \ t ts — t 



= H , (35) 



which satisfies (|21|l . 

We should note that any kind of finite large quantity could be chosen as the infrared cutoff. For instance. 
La could depend on dLp/dt and/or dLf/dt. The following model may be considered: 



1 / dLp dLf 



La Lp + Lf \ dt dt 

with the assumption that 



da dt 



Ha^ JO 



(36) 



(37) 
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is finite. Since c is a constant and Lp + L f = ac, from l|36|l . we find that H21|l is trivially satisfied. Therefore, 
in the model H36|l . as long as c in H37() is finite, an arbitrary FRW metric is a solution of H21() . 

The next step is to consider the combination of the holographic dark energy and other matter with the 
following EOS 

P=-P- fip) ■ (38) 

The holographic dark energy is given by (|20|) with L\ = Lp, Lf in or Lf in H28|l or H3U|) . Then from 
the first FRW equation 

First, we consider the case when f{p) =^ —(1 + w)p with constant w, that is, 

p = wp . (40) 

Then, by using the conservation of the energy 

0^^ + 3H{p + p) , (41) 

it follows 

p = poa-'^'+^^ . (42) 
By combining 139|l with (|42|l . one finds a solution 

{(2^^3(1+™) when w > —1 

2 
flo [ts - t) 3(1+") when w < -I 

3(i+„) 3 f 2(1 -c) , 1 f 2(1 + c) 1 

=;?l3(TTH''7l^aT^"7 ■ ^''^ 

c is assumed to be non-negative. Hence, in order for po^o^^^'"'' l|43|l to be positive, the constraint for the 
solution (|43|l appears: 

12 12 

Then, if w < —1, it follows c < 1. 

For another EOS (see also the previous section) 

f{p) = /op" , (45) 

the conservation of the energy (|41|l gives 



Defining a new variable q by 



smce 



p= 3/o(l-a)ln— ^ . (46) 
ao 



gEEln-, (47) 
ao 



!L^aH^ = H^, (48) 
at da aq 
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Ea. (|39|l can be rewritten as 



Thus, Eq.lgni) gives 




c 



-9 



H^~^P) f-^. (49) 



p = piq^, pi = {3/o(l-a)}— , (3=- . (50) 

1 — a 

As we are interested in the case when a, and therefore q, is large, we assume that H behaves as H ^ hiq^ , 
when q is large. From (|49|l . one gets 

(l-c')hl^^p,, (51) 

which requires < 1. At the same time, 

aoceV^-l / or eV ^-f J . (52) 

Hence, ifl — ^ <0or 1/2 < a < 1, a diverges a,t t — ts, which is similar to a Big Rip singularity 't^. 

Let us now consider the contribution of quantum corrections. Quantum corrections are important in 
the early universe when i ~ and curvatures are large. However, the quantum corrections are important 
near the Big Rip singularity, where the curvatures are large around the singularity at t = tg. Since quantum 
corrections usually contain the powers of the curvature or higher derivative terms, such correction terms play 
important role near the singularity. It is interesting to take into account the back reaction of the quantum 
effects near the singularity^^]. Note that holographic dark energy could be caused by non-perturbative 
quantum gravity. However, it is not easy to determine its explicit form. As it has been shown above, the 
Big Rip may occur also in a holographic dark energy model. On the other hand, the contribution of the 
conformal anomaly as a back reaction near the singularity could be determined in a rigorous way. In the 
situation when the contribution from the matter fields dominates and one can neglect the quantum gravity 
contribution, quantum effects could be responsible for the change of the FRW dynamics non-perturbatively. 
Hence, as in other dark energy models 6, 8] the Big Rip singularity could be indeed moderated or even 
prevented. 

The conformal anomaly Ta has the following form: 

TA^b(^F+ ^V^i?^ + b'G + b"V^R , (53) 

where F is the square of a 4d Weyl tensor and G is a Gauss-Bonnet curvature invariant. In general, with N 
scalar, iVi/2 spinor, iVi vector fields, N2 (= or 1) gravitons and iVnD higher derivative conformal scalars, 
the coefficients b and b' are given by 

N + 6N1/2 + l2Ni + 6II7V2 - 8iVHD 
120(47r)2 ' 

+ 117Vi/2 + 62iVi + 14117V2 - 28iVHD 

b = : — . (54) 

360(47r)2 ^ ' 

b > and 6' < for the usual matter except for higher derivative conformal scalars. We should note that b" 
can be shifted by a finite renormalization of the local counterterm i?^, so b" can be arbitrary. 

In terms of the corresponding energy density pA and the pressure pA, Ta is given by Ta = —pA + -^PA- 
Using the energy conservation law in the FRW universe: 



PA + 3iJ {pA + pa) = , 



(55) 
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one may delete pA as 

Ta = 'ApA - pa/H . (56) 

This gives the fohowing expression for pA'- 



PA = J / dta^HTA 



I to 

1 



J to 

(46 + 66") (h + 7HH + iH^ + UH'^h) ] . (57) 



The next step is to consider the FRW equation 



=PA + PA . (58) 



K 



The natural assumption is that the scale factor a behaves as H25(l when i '--^ or as ()27(l when t ^ ts- Then 
and pA behave as t"^ or {tg — 1)~^ but pA behaves in a more singular way as or {ts — i)~^- Hence, 
the power law behavior cannot be a solution of (|58|l . Instead, the spacetime approaches the deSitter space 
as m if^ . From ^ one finds 

When c < 1, which corresponds to the singularity (|25|l of the early universe, Ea. H59|l has solutions 



Thus, the singularity in the early universe may be stopped by quantum effects. The whole universe could 
be generated by the quantum effects. On the other hand, if c > 1 the only consistent solution of H59|) is 

H=\=0. (61) 

This shows that the universe becomes flat. Therefore due to the quantum matter back-reaction, the Big Rip 
singularity could be avoided. Hence, unlike to the model of previous section, the effective phantom phase 
for holographic dark energy leads to Big Rip which may be moderated (or prevented) by other effects. 



IV. DISCUSSION 



Summing up, we have considered in this paper two dark energy models both of which contain an effective 
phantom phase of the universe evolution without the actual need to introduce a scalar phantom field. Both 
models, which can be termed as decaying vacuum cosmology and (generalized) holographic dark energy, 
have a similar origin related with quantum considerations. However, the details of the two models are quite 
different. For instance, the decaying vacuum cosmology is caused by vacuum fluctuations and, under a very 
reasonable assumption about the decaying law, the effective phantom phase is transient, no Big Rip occurs. 
Moreover, the observable (effective) EOS parameter being the phantom-like one to comply with observational 
data is different from the real dark energy EOS parameter, which may be bigger than —1. At the same time, 
holographic dark energy is motivated by AdS/CFT-like holographic considerations (emergence of an infrared 
cutoff). Even without the phantom field, the effective phantom phase there leads to a Big Rip singularity, 
which most probably can be avoided by taking into account quantum and quantum gravity effects. 
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It looks quite promising that one can add to the hst of existing dark energy models with phantom-like 
EOS (some of them have been mentioned in the introduction) two additional models which exhibit the 
interesting new property that there is no need to introduce exotic matter explicitly (as this last is known to 
violate the basic energy conditions). Definitely, one may work out these theories in more detail. However, as 
usually happens, it is most probable that the truth lies in between the models at hand, and that a realistic 
dark energy ought to be constructed as some synthesis of the existing approaches which seem to possess 
reasonable properties. New astrophysical data are needed for the resolution of this important challenge of 
the XXI century. 
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